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We derive and apply an optical Bloch equation (OBE) model for describing collisions of ground 
, and excited laser cooled alkali atoms in the presence of near- resonant light. Typically these collisions 

lead to loss of atoms from traps. We compare the results obtained with a quantum mechanical 
complex potential treatment, semiclassical Landau-Zener models with decay, and a quantum time- 
dependent Monte-Carlo wave packet (MCWP) calculation. We formulate the OBE method in both 
adiabatic and diabatic representations. We calculate the laser intensity dependence of collision 
probabilities and find that the adiabatic OBE results agree quantitatively with those of the MCWP 
calculation, and qualitatively with the semiclassical Landau-Zener model with delayed decay, but 
that the complex potential method or the traditional Landau-Zener model fail in the saturation 
! limit. 

O ■ 32.80. Pj, 42.50.Vk, 42.50.Lc 

CI*; I. INTRODUCTION 

Collision dynamics of cold atoms in laser traps have been extensively investigated over the past few years. When 
the red detuning A from the atomic resonance frequency is large compared to the natural decay rate 7, we find a 
photoassociation spectrum of isolated bound vibrational-rotational levels in the attractive excited state potentials. 
This is now a highly developed subject and is fairly well understood, see the review Q]. In contrast, when A is small, 
on the order of 7, the mechanisms and rate coefficients of trap loss processes which result from photoexcitation of the 
. diatomic quasimolecule at long range are still rather poorly understood theoretically [^,^, in spite of the numerous 
J — . ' experimental studies of this subject Q. 

, The reason for this is twofold: real hyperfine structure introduces much complexity into the collision dynamics, 
' and the prominent role of excited state spontaneous decay during the very long time scale of the collision is difficult 
O i to calculate quantum mechanically. The number of degrees of freedom associated with the spontaneous emission is, 
' ^ . of course, infinite. Adiabatically eliminating these degrees of freedom leads to a mixed state representation that can 
not be described in terms of wavcfunctions but requires solving the Liouvillc-von Neumann equation for the quantum 
mechanical density matrix /3(_R, i?'; i) 
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R'- 1) = -^ [H{R)p{R, R'; t) - p{R, R'; t)H{R')] + Vp{R, R'; t), (1) 
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$-H where H{R) — T{R) + V{R) is the system Hamiltonian for kinetic energy T{R) and interaction potential V{R), and 
. 5^ , r is the decay tensor. Thus, the theoretical treatment of cold atom collisions serves as both prototype and paradigm 
for new constructs to treat non-equilibrium open systems coupled to reservoirs. Since the direct solution of Eq. (|^) for 
cold collision situations is beyond currently available computational resources [^|j7| , approximate methods for treating 
the collision dynamics in light fields must be developed. 

The methods currently available are the semiclassical local equilibrium model of Gallagher and Pritchard j|] or 
Julienne and Vigue |9|, the semiclassical dynamical Landau-Zener model s |10|-p^ , the semiclassical optical Bloch 
equation (OBE) method [Q, the quantum complex potential method |l^Jl5^0, and the Monte Carlo wave packet 
method of simulating the full quantum density matrix 1 11 - 13| . Although the latter is capable in principle of treating 



the full quantum dissipative dynamics for an arbitrary number of coupled states in arbitrarily strong laser fields, the 
method is extremely computer intensive, and therefore slow and impractical. The complex potential method can treat 
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fast many coupled channels quantum mechanically, including bound state resonances, but only in the limit of very 
weak laser fields where no more than one excitation and decay event per collision occurs. 

The semiclassical methods are very appealing because of their computational tractability, simple interpretation 
and physical picture of the collision. However, several fully quantum calculations |TT|-[T^ have shown that both the 
local equilibrium and semiclassical OBE methods (in the formulation given in Rcf. ||14|) give incorrect results by an 
order of magnitude or more for detunings of a few 7 or less for temperature T < 1 mK, depending on species. So 
far no practical theory exists for ultracold collisions for realistic atoms in a light field with A ~ 7 which is fully 
quantum mechanical and also capable of treating dissipation and decay. Therefore, there is not yet any satisfactory 
description of trap loss rates in the small detuning limit. For large detunings, A ^ 7, collision in a light field goes to 
photoassociation spectroscopy, in which isolated molecular bound vibrational levels are excited. Resonant scattering 
theory does then an excellent job of explaining the excitation rate p|,p^. 

Quantum calculations have shown that semiclassical methods may still be useful in characterizing cold collisions 
in a light field with A ~ 7 [pT|-p^. The local equilibrium model for cold collisions place a prominent emphasis 
on off-resonant quasimolecular excitation outside the region around the Condon point Rc- In contrast to them, a 
semiclassical picture based on localized Landau-Zener excitation near Rc with subsequent semiclassical evolution with 
decay inside Rc gives an excellent representation of the quantum dynamics for T near 1 mK, when compared with 
results obtained from quantum mechanical calculations. The Landau-Zener model only begins to fail near T = 1 fiK 
and at large intensities Therefore, since there still seems to be good opportunities for semiclassical models, 

we have revisited the OBE method, and provided a rigorous derivation of the velocity-corrected semiclassical OBE 
equations starting with the appropriate quantum mechanical equations (replacing the unsatisfactory approach used 
in Ref. We show that an adiabatic rather than a diabatic formulation of the semiclassical OBE equations gives 

quite good agreement with the quantum methods, even at saturation, in contrast to the poor agreement provided by 
the diabatic treatment. 

We find that the adiabatic OBE calculations are in good agreement with the MCWP results and time-independent 
complex potential results (used for weak laser fields where this method is appropriate), down to low collision tem- 
peratures. Only upon increasin g th e laser detuning above the onset of resonances due to bound state structure does 
the adiabatic OBE method fail pO[|. We also use an improved Landau-Zener model with dissipation jl^ which works 
even in the presence of strong saturation, where the complex potential method fails. This model offers a qualitative 
understanding of the strong field processes. The numerical comparisons are for the standard two-state model on which 
the quantum and other semiclassical models have heretofore been tested. Although these test calculations ignore the 
complex multistate structure introduced by molecular hyperfine structure, the hope is that semiclassical methods can 
yet be developed that are capable of treating the complexity of multistate collision dynamics in the presence of decay. 

This paper is constructed as follows. Section |l| presents the model of trap loss processes we shall use in order 
to test the methods developed and employed. Section |l| contains the derivation of the OBE method using the 



two different bases. Section IV describes the Monte-Carlo method which serves as the standard against which the 
approximate methods we use are compared. Section ^ describes the complex potential i neth od. Section VI develops 
the generalized Landau-Zener approach to strong laser field cold a tom collisions. Section VII presents the comparison 
of the numerical results from the various methods, and Sec. VIII contains a summary and conclusion. 



II. THE QUASIMOLECULE MODEL FOR COLD COLLISIONS 

The basic loss processes for atom traps due to cold collisions are fine-structure change (FS) and radiative escape 
(RE) these are demonstrated in Fig. |l|. Due to the low temperatures we can consider the collision of two atoms 

as internal dynamics of a diatomic quasimolecule. The simplest prototype model for the FS mechanism consists of 
three collision channels, i.e., quasimolecule states. The model described here is identical to that used in Ref. We 
ignore any rotational structure, i.e., only the s-wave is considered, but the model can be extended to higher partial 
waves. Here the three channels are: the ground ^5*1/2 + ^5*1/2 state channel labelled g, the excited + ^^3/2 

state channel labelled e and a probe channel (correlating asymptotically with ^51/2 + ^-Pi/2 state) labelled p. 

In the FS mechanism the system starts on channel g, and is later excited at the Condon point Rc to the channel 
e, which has an attractive potential. When the atom reaches the crossing between the potentials for the e and p 
channels, it may enter the p channel and eventually come out of the collision having gained as kinetic energy the 
energy difference between the ^S'1/2 + ^^3/2 and ^5*1/2 + states. This gain is large enough to propel atoms 

from the shallow trap. 

In the RE mechanism the quasimolecule decays back from the e channel to the g channel via spontaneous emission. 
If this decay does not take place too early, the acceleration on channel e will give the atoms enough kinetic energy to 
escape from the trap. 
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FIG. 1. The basic trap loss mechanisms. The figures show the quasimolecule potentials as functions of the interrmclear 
separation R. The corresponding asymptotic combinations of the atomic states are also given. In (a) we demonstrate the 
fine-structure change (FS) loss mechanism. The quasimolecule is excited from the ground state g to the state e, then it moves 
towards small R, where it is transferred to the p state at i? = Rin- Finally the atoms exit the collision on this state, sharing a 
kinetic energy increase equal to AE. In (b) the excited state e decays back to the ground state (g) after the atoms have gained 
enough kinetic energy {AE) to escape from the trap. This is the radiative escape mechanism (RE). If decay takes place too 
soon, i.e., at R > Rin, then the escape turns into heating due to insufficient increase in the kinetic energy (the trap depth is 
not exceeded). 
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For both mechanisms we need to find out the probability for the quasimolecule to reach a certain internuclear 
distance i?in while still remaining on the channel e. This can be obtained from calculations by monitoring the 
quantum flux Je{R) on channel e directly as is done when using the OBE method, the MCWP method and the 
Landau-Zener approaches (two-state case). Alternatively we can monitor the population transferred to the probe 
channel p, as happens in the complex potential method (three-state case) . 

We need a treatment that contains both the laser-induced excitation at Rc from g to e, and survival on e. Fur- 
thermore, for strong fields a proper model must allow any decayed population to be excited back to e if the decay has 
taken place in the vicinity of Rc- Only the MCWP method and the OBE approach can handle this reexcitation (also 
called population recycling) quantitatively. Further discussion on the subtle aspects of these loss mechanisms can be 
found e.g. in Refs. |^|J|. 

It should be pointed out that the OBE method, the MCWP method (as we apply it in this paper), and the Landau- 
Zener approaches are "one-way" studies. We only consider the flux going in, but do not allow for any outcoming 
flux. For estimating the FS and RE loss this is adequate as long as the detuning of the laser is about 1-10 atomic 
linewidths. Then the role of bound states in the loss mechanisms is not too important yet; see Ref. for more 
discussion. For strong fields the power broadening also diminishes the role of the bound states; most of the loss is 
due to processes associated with the first passage of the critical point -Kin on the channel e. 

We have earlier in a short article (is) presented the results obtained with the adiabatic OBE method and the MCWP 
method in the strong field regime, and discussed the physical implications of the results (the lack of saturation in trap 
loss when excitation becomes saturated). In this paper we study in detail the various theoretical approaches, present 



the general derivation of the OBE equations and their application to the two-state case (Sec. Ill) and show how to 
extend the Landau-Zener approach to the strong field regime as suggested in Ref. in the case of optical shielding. 
The model Hamiltonian is 

u-i[K)L+y hA + UeiR) + Uc{R,l) J ' ^' 

where T(R) is the radial kinetic energy operator, Ug{R) is the ground electronic state which behaves asymptotically as 
Ug{R) = C6/i^^ Uc{R, I) = h^l{l + l)/2^iR^ is the centrifugal potential, fjL is the reduced mass of the quasimolecule (we 
assume a Cs2 system), Ue{R) is the excited state potential correlating asymptotically to ^5*1/2 -I- ^P3/2 alkali atoms 
and behaving asymptotically as —C^/R'^, and A is the detuning from resonance, HA — E{P3/2) — E{Si/2) — tT-LO- 
Here to is the laser frequency, and the laser-induced coupling is described by the Rabi frequency $7. The g and e 
channel potentials cross at the Condon point Rc{A) where hA — Ue{Rc) — Ug{Rc)- This crossing occurs at large 
internuclear distances. The values used for potential parameters are C3 — 20.30 e^a^ and Cq = 6.40 x 10^ e^a§, where 
e is the electron charge and oq — 0.529 A is the Bohr radius. Note that we express all energy and angular frequency 
parameters in frequency units. 

In studies using different laser parameters we have selected for a suitable inner distance i?in = 143 ag, although the 
methods that we use allow the determination of Je{R) for all values of R. The atomic excited state has a linewidth 
7at = 5.13 MHz, and we have taken the molecular linewidth to be 7 = (4/3)7at, independent of R. And as mentioned 
before, we consider only the case ^ = 0. 

A reasonably complete description of cold atom collision dynamics in laser traps can be obtained via the time- 
dependent density matrix p(R,R';t) satisfying the Liouville equation ([^). The probability of reaching the inner 
region on the excited state potential (or the particle flux in the inner region on the excited state potential) can be 
determined directly from the diagonal density matrix elements pu^Rim R'm] t). However, the direct numerical solution 
for the density matrix is presently beyond our capabilities for the cold atom collision problem. Instead we shall use 
various approximations to solve for the dynamics. The time-dependent MCWP approach basically includes all the 
physics contained in the Liouville equation, and it will provide the standard against which all other methods are to 
be judged. 

Our model is a simplified representation of the true collision situation. Experiments have shown that inclusion of 
hyperfine structure is necessary to properly treat the collisions of laser cooled alkali species |^,|l^ . It is exactly for this 
reason that it is so important to develop simple and accurate approximate numerical methods that can conceivably 
be used on problems including hyperfine dynamics where a large number of channels is required to treat the manifold 
of the hyperfine states. 



III. OPTICAL BLOCK EQUATIONS 
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A. General situation 



In this section we derive the semiclassical optical Bloch equations. We start with the LiouviUe equation (||) and use 
Wigner function description into which we introduce semiclassical approximations. If we write Eq. (|^) by components 
we get 



^[■K,fc(i?)Pfcj(i?, i?'; t) ~ pM i?'; t)Vkj{R')] + ihY^T,jkiPki{R, R'; t). (3) 



kl 



Here we have explicitly written out the kinetic and potential energy parts of the Hamiltonian H; the term V{R) 
contains the potentials for the internal states of the quasimolecule and couplings between them. We assume that V 
has no time-dependence as we have eliminated the oscillating laser field terms using the rotating wave approximation 
and an appropriate phase shift (here V contains the potentials U and couplings h,^l). 

It should be pointed out that our description is time-dependent, so instead of the boundary conditions used in the 
time-independent scattering theory we have an initial value problem, i.e., we solve Eq. starting at t — to with 
some initial density matrix p{R, R'; io)- If Eq. (^) corresponds to a closed system (no decay out of the selected set of 
levels), then p can have steady state solutions. 

Typically one takes as the initial state the steady state result corresponding to atoms being well separated, with 
the quasimolecule potentials being flat over the distance that the system moves within the time it takes to establish 
the steady state. In other words, the molecular potentials do not impose any dynamics that would interfere with the 
steady state formation. 

In time-independent scattering theory the initial conditions can not be stated in terms of diabatic states if couplings 
between the states do not disappear asymptotically, as is the case with laser-induced quasimolecule processes. However, 
since the system at large R evolves quickly into the local steady state, which is independent of the selected basis states, 
there is no a priori reason to regard the adiabatic basis better than the diabatic basis. In practice one tends to choose 
the diabatic basis, because it allows a simple description of the spontaneous emission processes. Furthermore, as 
discussed later, we can select any initial state in any basis, if we allow the system initially enough time to reach locally 
at large R a steady state before the spatial dependence of the quasimolecule potentials will couple the steady state 
formation and molecular dynamics. 

We assume for simplicity that in our current description F is independent of position, but it is easy to extend our 
treatment to allow R dependence in F; such dependence can easily arise if retardation effects are properly included 
to the quasimolecule potentials and lifetimes. We redefine our spatial coordinate system by writing R = r + q/2, 
R' = r — q/2, which transforms the kinetic term: 

0^ 



9i?2 dR'^ drdq 

The density matrix p(i?, i?'; t) contains information about the spatial coherences in the system. In order to calculate 
quantum fluxes at some interatomic distance R we do not need all that information, but only the spatially diagonal 
elements p{R, R; t). However, the evolution of these diagonal elements depends on the off-diagonal p{R, R'; t) elements. 
By using the Wigner function 

f°° 11 
Wij (p, = y dQ exp{-ipq/n)p,j {r + -q,r - -q-t) (5) 

we can include the spatial coherences and yet effectively work with the spatially diagonal terms only. 

First we apply the Fourier transform given in Eq. ^ on both sides of Eq. (|^) in order to obtain the equation of 
motion for the Wigner function: 

dqexp{-ipq/h) \ ^ J2^V^kir + q/2)pkj{r + q/2, r ~ q/2; t) - p,k{r + q/2, r - q/2; t)Vk,{r - q/2)] 
r'^ fc 

Y,^^3kiPkl{r + q/2,r ~ q/2;t)\ . (6) 
kl J 
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One should note that above we have appUed integration by parts in order to replace — with p; this requires that 

\m\q^±oo p{r + ^q, r — ^q; t) = 0, i.e., that the spatial coherences disappear as we move away from the diagonal — this 
is a reasonable assumption. 

If we integrate the Wigner function over momentum p we get the spatial probability distribution, which we can 
define as 

~P-o{r,t) = —j dpW,jip,r;t). (7) 

The quantity pij{r, t) equals Pijir, r; t), as can readily be seen by substituting the expression for Wij(v, r; t) in Eq. ^ 
into the right hand side of Eq. (0) and carrying out the integration over p. So, by integrating Eq. (g) over p we get 
the equation of motion for pij{r,t): 



9 ~ , 15 



1 f°° 1 1 

— / dppW^j(j),r;t) = —^[V,kir)pkjir;t) - p,k{r;t)Vk-i{r)]+^T,jkiPki{r,t). (8) 



Here we have used the fact that there is no p-dependence on the right-hand side of Eq. (|6|) , so that if we perform the 
momentum integration first, we obtain a 5{q) function, and thus the integration over q merely sets q = 0. 

If a convenient method of evaluating the kinetic term were available, we could use the result (H) to obtain exactly 
the diagonal elements of p at given r and t, pii{r,t), which is the probability of being in channel i at position r and 
time t. However, since we do not have any exact methods for calculating the second term in Eq. (||), we shall estimate 
it using a WKB approach. The WKB approximation for the density matrix element is given by 

ai(R;t)a*(R';t) 

p,,{R, R'; t) ^ ) \''\^' eMiim) " PARW^ (9) 

y^p,{R)pj{R') 

where Pi{R) is the local classical momentum in state i, (3i{R) is the action at R, 

(3^{R) = I dxp^ix), (10) 



and ai{R,t) is the amplitude factor for the WKB wave. We consider only the incoming wave, which fixes the sign of 
the P terms, and assume that and pi do not depend on position very strongly. Furthermore, we assume that the 
classical momenta Pi{R) are non-zero and real. It should be noted that by introducing the classical momenta we have 
made our equations energy dependent as well, since 



V^{r) = ^2p{E~[Vu{r)^Vu{^)]}, (11) 

where E is the asymptotic energy for the WKB wave (equal to the asymptotic relative kinetic energy of the colliding 
atoms). 

It should be noted that our definition (|l^) of the classical momenta is clearly different from the one encountered in 
the traditional scattering theory, if we consider the asymptotic situation. In the time-independent theory the channels 
(states) are typically either open or closed, depending on the collision energy, i.e., their classical (WKB) momenta are 
asymptotically either real or imaginary. This is because they are defined as p^^^^^ (r) = yj2p{E — [Vii(r) — Poo(oo)]}, 
where i = corresponds to the channel of the ingoing wave, defined by the asymptotic boundary conditions. 

The difference here is due to the presence of the relaxation terms in Eq. (|^), and is required by the asymptotic 
situation. In the time-dependent treatment we have Pi(oo) = ^2pE, which is independent of the state label i. This is 
because asymptotically we have a steady state formation which is not coupled to the dynamics because the potentials 
are flat. For simplicity we base our following discussion on a two-state system. Assuming that the excited state and 
ground state populations had different asymptotic momenta, the steady state formation (the cycles of excitation and 
decay) quickly mixes these populations and eventually the distribution of momentum on the ground state and the 
excited state would be exactly equal. In other words, because of decay we can have asymptotic population even on 
a closed channel, but this population is a steady state reflection of the ground state population and must have the 
same classical momentum. Although in the above discussion we have assumed the diabatic basis, it is quite valid in 
the adiabatic basis, where the relaxation leads to a similar asymptotic mixing between the channels. Fin ally, a s we 
fix the asymptotic situation by using Eq. (^), we introduce other problems, which will be discussed in Sec. IIIC. 



Next we insert Eq. (||) into Eq. (|5|), and use the result in Eq. (^). The exponential part of the integrand can be 
expanded around r: 
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exp{z[A(r + iq) - f3j{r - ^q)]/h} ~ exp{z[A:(r) - P,ir)]/h + t^[p,{r) + pj{r)]q/h + 0{r^)}. (12) 
Then we apply the stationary phase method to obtain 

1 r°° 1 

— J dppWij{p, r; t) = - [pi{r) + Pjir)]pij {r, t). (13) 
Substitution of this expression into Eq. (^ yields 

—Pij{r,t) + ^[p^{r) + Pj{r)]—p.ij{r,t) = — Y^[V^k(r)pk3{r;t) - p,k{r;t)Vkj{r)] +'^T^jkiPki{r,t). (14) 

^ k kl 

Here we have assumed that the classical momenta vary so little with r that they can be taken outside the derivative 
term. 

Our aim is to find the total incoming quantum flux at each position r (integrated over all times), and thus we are 
not interested in the actual time dependence. This simplifies our model to a great extent. Now the total fiux can be 
obtained as a steady state result from Eq. (p^. Since our model corresponds to a "one-way" situation, the steady 
state flux at r is equal to the total flux that has passed that point. In the steady state the time derivative in Eq. (|lj) 
vanishes, and we can replace py (r, t) by {pki)ss{'f)- Furthermore, it is convenient to define the quantity 



f^kiir) = \/ Pk{r)pi{r){pki) ss{r) (15) 

whose diagonal elements give the flux in the various states at position r. We shall call this quantity the semiclassical 
density matrix. The equation of motion for it, within the validity range of the WKB approximation, is 



where 



^ kl 



1 



Pt{r)pj{r)^^^'^ 



Pk{r)pi{r) 



[Lijki +Yvjki\<Jki{r), (16) 



Ujki = -^[V,k{r)6fl - 5,kVij{r)]. (17) 



This completes the general derivation of the semiclassical optical Bloch equations using the Wigner distribution. The 
set of equations ( |l6| ) were obtained by making a semiclassical approximation, which focuses on classical paths by 
virtue of using the WKB approximation. Only when the semiclassical approximation is valid (when the de Broglie 
wavelength is smaller than the region where the potentials are varying) will this approximation be meaningful. 

The above derivation was done in the diabatic representation of the quasimolecule potentials. The decay term F has 
a simple form in this representation, where the electronic states are also independent of position, and the internal states 
are directly coupled by the standard dipole term. In the adiabatic representation the electronic wavefunctions vary 
with the internuclear coordinate R. In this presentation the laser-induced couplings and the decay term become clearly 
R dependent. We can move from the simple diabatic representation into the adiabatic one, in which the potential 
matrix V (which contains the radiative coupling) is diagonal. The transformation matrix C{R) is R dependent, 

^ Cr-^\R)Vki{R)Ci,{R) = 6,,E,{R), (18) 

kl 

where E{R) is the diagonal eigenvalue matrix that gives the field-dressed quasimolecule potentials. We can now write 
the diabatic semiclassical density matrix in terms of the adiabatic semiclassical density matrix pf^iR, R'] t): 

p,,{R,R'-t)^Y.^,^^^)pli{R,R'-t)Ci^{R'). (19) 

kl 

Next we insert this pik{R, R' , t) into Eq. (^), and proceed as in the diabatic case. When making the WKB approxima- 
tion and using the stationary phase approach we assume that the matrix elements Cij (i?) are slowly varying functions 
in position. Thus, by sandwiching the whole equation between C^^ and C we obtain eventually 



^[rf(r) +p^(r)]^at^ir) = ^[E,ir) E,{r)K, ^ [Q^,kl^[ptir) + pfir)] - ^ki] <^kiir) 



Piir)pUr] 



-.1/2 



Pl{r)Pi{r) 



(20) 
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where the classical momenta p° (r) are now defined using the adiabatic potentials Ei (r) in Eq. (jT^) . The decay term 
transforms as 



-f-i 



klmn 

The non-adiabatic couphng that arose when we evaluated d/dr[Cp°-{r, t)C^^] is given by 

Qtjki = ^ 

= E 

771 

= Ak{r)Sji - S,kAij{r), 



C-J:{r)^ank(r)S,i - S,kCr^{r)-^ani{r) 



(21) 



(22) 



where A,fe(r) - E™C'r„i(r)fC„fe(r). 

Another method for deriving the adiabatic OBE equations of motion involves using the half-collision matrix 
method |Q. This method yields the same result as given by Eqs. ( po|) and (^) for the Hamiltonian part of the 
dynamics, but can not be used to derive the decay part of the adiabatic OBE equations, since the half-collision 
method does not incorporate the dissipative dynamics due to spontaneous emission contained in the density matrix 
treatment. 



B. The two-state case 



We assume that our quasimolecule has only two states, one ground state (1) and one excited state (2), with 
potentials Vn and ¥22- In our trap loss model these states are as shown in Fig. |[ The excited state has a constant 
width 7, and the off-diagonal density matrix elements (pi2 and p2i) have the width ^7. Thus we have 

ril22 ~ 1, r2222 = —J, r2121 = ^2^^' -'^1212 = ^2^' (23) 

and rest of the elements of F are zero. Using V12 as the coupling between the states, we get the diabatic equations 

~a2i)+7— (24) 

(25) 

) +»7 « 2Vi2 f \ /ORA 

Cri2 - ■ 1/ — 0-22 - W — CTii (26) 

2 n Vi + V2 \\j V2 \j Vi ) 

(27) 

Here we have used the classical velocities Vi[r) — piir)/ p,. 

In order to move into the adiabatic frame we need the transformation matrix elements, Cij{r). In a two-state 
system it is convenient to define 

^=itan--^. (28) 

Z V22 — Vii 

The transformation matrix is given by 

C{r) = ( ^" = ( ''"[^i ) , (29) 

^ ^ \C2i C22 J \- sin{9) cos{e) y ' ^ ' 

and the inverse transformation is obtained from the relation C^^{9) — C{^6). Thus we get 

dr~ dr 86^ dr\-l Q)- ^^^"^ 



da 11 


I V12 , 


dr 


" % 1 ( 

ft V^it;2 


d(T22 


dfTii 


dr 


dr 


dai2 


I 2(T/ii - 


dr 




d(T2l 


_ dal2 


dr 


dr 
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FIG. 2. The two-state model for trap loss collisions. The solid lines describe the bare (diabatic) quasimolecule potentials 
near the Condon point Rc- Here A = 7, so Rc — 2930 ao. The dotted lines describe the field-dressed (adiabatic) potentials. 



For convenience, we define D = This gives us 

Q1121 = Q1222 = A12 — D, Q1112 ~ = D, Q1211 = —A12 = —D. 

After a little algebra, we get the semiclassical optical Bloch equations in the adiabatic representation: 



dot, 
dr 

dr 
da^ 
dr 



da^ 
dr 



-D- 



Ki+f^?2)-7 



'11 



^22 



/^a I a \ 




D ^ ^ ^ g° 



sc(l 




(31) 

(32) 

(33) 
(34) 



Here we use the notation s — sin(6') and c = coa{6). The velocity factors and wf are defined as in Eq. (jTl 
El and i?2 as the appropriate potentials. At this point we relabel r with R. 



(35) 
using 



C. The implementation of the two-state case 



The optical Bl och e quations (p^-(p7|) and (p^-(|35|) can be solved numerically using various methods. As will be 
discussed in Sec. VIl, the diabatic formulation fails when T < 1 mK for all coupling strengths $7. If we use the 
adiabatic equations (|32|)-(|3^), with the velocity factors given by Eq. (11), we find a good agreement with the MCWP 
results for all temperatures at large couplings, but an increasing deviation with decreasing f2. Of course, we would 
like to have an approach which is good for all T and fi. Thus we need to look into the issue of velocity factors in 
detail. 
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The velocity factor V2 for the upper adiabatic potential (see Fig. ^ is the source of the problem here. When the 
atoms approach each other, a steady state is quickly formed and maintained until the system gets very close to Rc- 
As discussed earlier, this gives us the advantage of a basis-independent initial state. At this point the population 
on the upper adiabatic state 2 is a steady state reflection of the population of the lower adiabatic state 1. Thus 
when the system reaches the Condon point the probability flux on both adiabatic states ought in fact have the same 
momentum. Thus, if we use the velocity factor given by Eq. (11) for the upper adiabatic state, we get the right velocity 
asymptotically, but the wrong velocity at the crossing. We can correct this by redefining V2{R) so that tif = y/2E/^ 
for all R. 

In other words, because the OBE approach mixes concepts from time-dependent theory (steady state formation) 
and time- independent theory (WKB wave functions), we must sometimes improve its performance by such a tun- 
ing. Here we need to give the WKB wave functions asymptotically the velocities which describe the steady state 
situation correctly, but at the crossing the dynamics dominates, implying that the classical momenta used in the 
time-independent theory would provide a more accurate description. Deviations are visible in the limit of the weak 
excitation, because then the excited state population is dominated by the small steady state contribution which, due 
to the large momentum given incorrectly by Eq. (11), has a good chance to survive on the excited state until R = Rm- 
This population overwhelms the contribution from the dynamical excitation, and gives thus incorrect results. How- 
ever, as will be shown in Sec. |VII| , we now have a method that can predict correctly the probability to reach small R 
on channel e for any practical Q, and for the temperatures in the cold collisions regime. It should be noted that an 
extension of our approach to ultracold temperature regime, i.e., below the recoil limit, is not likely to succeed, as the 
semiclassical viewpoint fails in this regime [^,^ . 

The previous formulation of optical Bloch equations was done in the time frame rather than in the position 
frame. The transition between these two frames was performed by introducing a reference trajectory vo{R), which 
mapped R to t. However, as noted in Ref. [14], this reference trajectory is not really needed as we can do the 
calculation in the R frame altogether. 

In practice when solving the optical Bloch equations we can set all aij{Ro) — 0, except selecting one state for which 
CiiiRo) — 1. The value of Rq is set suitably large for the atoms to be well apart and potentials flat. Since the classical 
trajectory couples R and t, evolution in R corresponds to evolution in t and we find that the system has evolved into a 
steady state distribution of the ground and excited state populations after moving a relatively short distance towards 
smaller R. This corresponds to numerically determining the asymptotic state state populations and coherences. 



IV. MONTE-CARLO SIMULATIONS 

A direct wave packet treatment of Eq. (|^) with numerical methods is possible [^,0. However, since one has to 
operate with a two-dimensional spatial grid, the memory sizes currently available in computers strongly limit the use 
of this approach: only models which are simplified in the extreme can be studied. In cold collisions the acceleration 
of the initially slow wave packet on the steep excited state potential surface forces us to use a large two-dimensional 
momentum space while at the same time good momentum resolution is needed to define adequately the narrow initial 
wave packet for low temperatures. Similar demands are set for the position space as well; for more detailed discussion 
see Refs. ]^Jll|. Hence we need many grid points in order to span properly the required regions in both the momentum 
and position spaces. 

It may be feasible to avoid some of the computational problems by using grid sizes and resolutions which are 
adaptive; one might utilize the rather deterministic behavior of the wave packet by altering the computational grid 
properties either as a function of position or time. Then the straightforward swapping between momentum and 
position representations using fast Fourier transforms is, however, usually lost. We have chosen to approach the 
problem from another angle. The Monte Carlo wave packet (MCWP) method allows us to treat Eq. ([|) numerically 
as a one-dimensional problem. Unlike other approximative methods this approach does not adapt any concepts from 
classical mechanics and, therefore, it is not a semiclassical tool but a fully quantum one. Hence it can be used as a 
benchmark for the different semiclassical methods described in this article. 

By using the MCWP method we can greatly diminish the limitations set by the available computer memory. 
However, this gain is partially reduced by the increase in the time required by the computation. Hence the MCWP 
simulations are quite time-consuming, and therefore the need for other approaches is quite acute. There are several 
Monte-Carlo approaches available currently, and we use the Dalibard-Castin-M0lmer version |^^, adapted to wave 
packet problems in the manner described in Refs. ||,|ll|,|l|. We shall give here a brief description of the method, but 
keep the main emphasis on aspects related to the particular system studied in this article. 

In the MCWP method one does not directly solve the time evolution of the density matrix itself. Instead, we look 
at the time evolution of the state vector 
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(36) 



where and ^'e are the ground and excited state probabiUty ampUtudes, respectively. In this model the spontaneous 
decay appears as random quantum jumps during the time evolution of the state vector. Hence each time we solve the 
time-dependent Schrodinger equation we obtain a unique state vector evolution, also called a wave packet history. One 
can form a finite ensemble of such histories and calculate ensemble averaged expectation values for physical quantities. 
These values are approximations to those provided by the full density matrix treatment. The accuracy of the ensemble 
averages tends to increase with the number of members in the ensemble, and in the limit of an infinite ensemble these 
averages and the density matrix results become equal, as shown e.g. in Ref. p3|] . So, we expect that by accumulating 
ensemble members we can eventually reach a suitable accuracy at some finite ensemble size. The accuracy to be 
expected of the method is decribed in Ref. . For wave packets in cold collision problems the appropriate ensemble 
size seems to be roughly 50 members, assuming that all the histories are very close to each other in the phase space 
for all times. This is quite true for the excited state survival studies related to attractive excited states. In general 
such localisation in the phase space for all times is necessary for the success of the semiclassical approaches. 

We start the calculation of the state vector evolution from some initial state, which in our case is a Gaussian 
wave packet on the ground state moving towards small R. The wave packet in general describes the probability to 
find the two colliding atoms at certain relative separation R, and its components |\E'g(i?, t)p and |^'e(-R, i)P contain 
the additional information how this probability is distributed between the ground and excited states. A Fourier 
transform of the state vector ^'(i?, i) takes the system into the momentum representation. We set the initial phase 
of in position representation such that the wave packet starts with a mean momentum (p) which corresponds to 
the temperature of the cloud of cooled and trapped atoms. The width of the wave packet is chosen so that it remains 
relatively narrow in both representations. We cannot, of course, violate the Heisenberg uncertainty relation, so it 
is impossible to have infinitely narrow packets in either representation. It should be pointed out that apart from 
satisfying the Heisenberg uncertainty relation the width of the wave packet is not related to any of the macroscopic 
quantities of the physical situation which we try to simulate. 

Initially the wave packet is located far from the crossing so that a steady state between the ground and excited 
states can form before the wave packet reaches the interaction region where the dipole-dipole interaction makes the 
laser resonant with the molecular transition. The time scale for the formation of the steady state is roughly 3-5 times 
the decay time scale l/7moi |Q, assuming that the local detuning does not change much over the distance covered by 
the wave packet during that time. As discussed before, the steady state formation allows us to put the initial wave 
packet on the diabatic ground state even when the laser-induced coupling is large asymptotically, because the final 
steady state is independent of the initial state. Indeed, we could even place the initial wave packet on the excited 
state, and yet the wave packet approaching the crossing would still be the steady state one. 

The state vector corresponding to the initial state is stepped forward in time with its evolution determined by the 
Schrodinger equation. Various numerical methods can be used, and we have applied the combination of split operator 
approach with Crank-Nicholson and Runge-Kutta algorithms, described in detail in Ref. [Q. In the MCWP method 
one uses an effective Hamiltonian, 



where H is the system Hamiltonian, 7 is the decay rate for the excited state population and (t+ and cr^ arc the 
standard spinor raising and lowering operators, respectively. 
For each time step t ^ t + St we calculate the jump probability 



where Pe{t) is the current excited state population. By rewriting 6s as dPe we would end up with the standard 
exponential decay exp(— 7t) of the excited state population. Now, we continue by comparing the jump probability Ss 
with a random munber rj Cz [0, 1]. A quantum jump occurs when 77 < 6s] this is usually the less likely situation since 
the basic assumption in the derivation of the MCWP method is that Ss ^ 1 all the time (guaranteed by choosing 
St <C 1/7). When a jump occurs one simply replaces '^g{R, t -\- St) with 5'e(-R, t + St), and then sets ^'e(i?, t + St) = 0. 
The occurrence of the jump corresponds to the observation of a fluorecence photon, which reduces the wave function: 
before the jump it had to be in the excited state, and after the jump it must be in the ground state. The important 
aspect is that as the jump takes place the position and momentum properties of the excited state component of the 
state vector are transfered to the ground state component. This is the source of radiative heating, among other things. 

Both the evolution under H^q and the quantum jumps reduce the norm of the state vector ^. Hence after each 
time step the state vector is renormalized to unity, even if a jump does not occur. It should be noted that if we had 



^^cff 



u -^T + - 
H — I — a a , 

2 



(37) 



Ss^^P,{t)St, 



(38) 
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Heff = then we would always observe a quantum jump for a system with non-zero Pe if we wait long enough. 
For cases where Pe < 1 this would be wrong, since there is a non-zero probability that the system never was on the 
excited state. In the weak field limit is always very small, so most of the ensemble members correspond to the 
time evolution under the non-Hermitian Hamiltonian with no interruption by jumps. Then the wave packet approach 
reduces to a time-dependent version of the complex potential approach. Therefore a single ensemble member becomes 
a reasonably accurate approximation to the density matrix result. We have used this property in our weak field 
study , and have thus verified that the diabatic formulation of the OBE method does not work properly at low 
temperatures, but the Landau-Zener approach and the complex potential method can be used instead. 

Although the MCWP method allows us to use relatively large grids, the strong change in the excited state potential 
corresponds at our probing distance i?in = 143 oq to kinetic energies which are beyond the numerical treatment. 
Basically it becomes impossible to correctly track the relevant quantum mechanical phase term exjp{—iESt/h), where 
E is the kinetic energy of the wave packet. Hence we must cut the excited state potential change by making it flat 
for i?in < R < Rcuu in our studies we have used the value i?cut = 512 ao- At i?cut we are basically left with the 
exponential decay of the excited state population because of the large local detuning. Hence we can take the wave 
packet result for Je(i?) from R — Rent to i?in by multiplying it with exp(— 7ttians), where itrans is the time it takes to 
go from i?cut to i?in along the classical path determined by the local velocity. In fact, the same approach is also applied 
when the OBE equations are solved numerically: otherwise the adequate determination of the term exp[—ip{R)6R/?i\ 
would require unpractically small values of the spatial grid spacing SR (here p{R) is the local momentum at R) when 

R < Rent- 



V. COMPLEX POTENTIAL CALCULATIONS 



In the complex potential method one adds a complex term on the excited state potential in order to describe decay 
out of these states [^0|jl^,|l^. This approach does not allow any reexcitation, and is thus not appropriate for strong 
field studies directly (by including the photon states explicitly one might improve the model although this would 
drastically increase the number of channels required to solve even the case of two quasimolecule states [p^). 

In this method one simply uses the Hamiltonian 

/ Ug{R) + UciR,i) nn \ 

H = T{R)+\ nn hA + UeiR) + UciR,i)~inj/2 nn^p , (39) 

V nn^p Ep + Up{r) + Uc.{R,l) J 

and solves the time-independent Schrodinger equation 

d'^F(R) 2u, , , , 

+ -^[E- U{R)]FiR) = 0, (40) 

where F is the three-component state vector for our model and E is the asymptotic collision energy. 

For this method we have explicitly included the probe channel p to our Hamiltonian: Up{R) is the corresponding 
potential. The probe channel p crosses the e channel potential at i?in. Because of the disparity between Rc{A) and 
Rin the outer zone excitation process is in practice well separated from the inner zone process. The coupling figp 
depends on the nature of the coupling of the excited state and the probe channel. 

We have used the invariant imbedding method ||l^,^ (in the diabatic representation) to solve the above close 
coupling equations in a form that directly computes the S matrix elements, Sgp and Sep, from which the quantum 
flux Je{R), i.e., the quantum mechanical probability of reaching the inner zone, is determined. In the complex potential 
calculations we use mainly the same values for parameters as in the two-state model, given in Sec. H. However, here 
the excited state potential, Ue{R), is taken as a numerical spline having a minimum energy of —182 GHz at i? = 72 oq 
and an asymptotic behavior of —Cz/R^- We take the probe potential to be Up{R) = C^/R^, with Cf — 7.260 c^Oq. 
All potentials have repulsive inner walls so the i? < region is non-classical. The other parameters are Ep = —3.0 
GHz and n^p = 1.0 MHz. 

As we are now working with the time-independent scattering theory, we need to consider boundary conditions 
instead of initial ones. Since the model potential contains non-vanishing off-diagonal elements, for strong fields there 
is a clear mixing of states at large R. Therefore the boundary condition for the complex potential approach has to be 
defined in terms of the field-dressed, i.e., adiabatic states. The transformation between the different bases then takes 



place as described in Sec. [II. We assume an incoming wave (corresponding to the asymptotic collision energy E) on 



the lower adiabatic channel formed by states g and e (state 1 in Fig. 
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VI. LANDAU-ZENER APPROACHES 



Since the inner and the outer crossings are isolated in distinctly different regions, we can write the loss probability, 
i.e., the probability to exit on channel p as 

Pp = \Sgv? = Pep(i?in) Je(i?in), (41) 

where the probability Pgp = I'S'epP measures the quantum probability of the inner zone e p process due to traversing 
the inner curve crossing once in both directions. Here Rin is the location of this curve crossing. When using the complex 
potential method we obtain the flux simply by using the relation Je(i?in) = |5'gpP/Pep(i?in)- The probability Pep is 
almost completely insensitive to the laser intensity / and to the collision energy E for small detunings. 

In the limit of large detuning and small laser intensity it is safe to assume that the excitation becomes localized to 
the Condon point Rc- In this limit the dynamical models and the local equilibrium model tends to agree. In local 
equilibrium models one assumes that the motion of the atoms is very slow compared to the steady state formation, 
and thus the steady state formation dominates This leads to a picture where off-resonant excitation is important. 
One can express the local excitation in a two-state model in the steady state regime as 

" A(P)2 + 2f]2 + (^/2)2 ' (42) 

where A(i?) is the local detuning, hA{R) — Ue{R) — Ug{R) -I- fiA. As a final stage in the local equilibrium model one 
weights the results with Eq. (^) and integrates over the position coordinate R. 

In the dynamical models it is assumed that as the system approaches Rq the motion and thus the change in the 
local detuning become fast compared to the steady state formation, and thus the excitation becomes a dynamical 
process which is localized to a region near Rc ■ The dynamical excitation can then be described with the Landau-Zener 
curve crossing model |10-1^]. In steady state models for large detunings the integration over the linewidth function 



( [42| ) becomes like a 5 function which singles out the Condon point, and thus the two viewpoints agree in this limit. 
The MCWP simulations have so far supported the dynamical picture over the local equilibrium picture; for further 
discussion see Ref. |]. 

In the weak field limit we can assume that the excitation and subsequent decay are uncoupled, and reexcitation is 
negligible. Then we can write, using the Landau-Zener model, the expression for the flux on channel e as 

J^'^^'iR) = SeiR, Rc;E, 7)Plz = S^iR, Rc;E, - exp(-27rA)]. (43) 

Here 

aVg(Rc;E) 

where Vg{R; E) is the classical velocity associated with the ground state at position R when the collision energy is E, 
and 



dUe{R) dUg{R) 



dR dR 



(45) 

R=Rc 



In other words, Plz is the one-way Landau-Zener probability of undergoing a transition from channel g to channel 
e at the Condon point Rc{A), and Se{R, Rc] E,'-f) is the survival probability from Rc to R. We can calculate the 
survival probability by assuming a classical trajectory combined with exponential (Weisskopf-Wigner) decay: 



' = exp(-7tei); Rc; E) = (46) 

Jrc Mr) 



where v^, (R) is the classical trajectory velocity for the excited state (in the diabatic formulation) . This is the Landau- 
Zener model with decay (LZD). 

The above model fails when excitation and decay do not decouple, which happens at strong flelds due to reexcitation 



of decayed population 1 1 1 |l^ . We can think of reexcitation as a process which delays the start of the exponential 
decay. Reexcitation takes place mainly within some region around the Condon point. We can deflne an interaction 
region for which A(i?) < fl. By making the simple assumption that exponential decay can take place only outside 
this region, we can rewrite the id in Eq. (E^) as 



13 



f dR 

where Rq is defined by the relation A(i?n) — Q; Rq < Re- The modified survival term depends now on the laser 
intensity / (/ (x fi^). We call this approach the Landau-Zener model with delayed decay (LZDD). Obviously the 
model can only give qualitative predictions, especially as the concept of the edge of the interaction region is not 
well defined. However, it gives a good intuitive undestanding why Je{R) for small R keeps increasing even when the 
excitation saturates to unity (and thus the LZD prediction saturates) . This picture agrees qualitatively with the 
results from the OBE and MCWP calculations, as will be shown in the next section. 



VII. COMPARISON OF METHODS 



Typically the excited state flux Je(R) shows oscillations at i? < Rc in the bare state picture. This is demonstrated 
by the MCWP results given in Fig. The oscillations are due to the coherences between the two states, established 
near the Condon point. As the coupling Q increases, the situation becomes increasingly adiabatic and the oscillations 
disappear. At the same time the asymptotic (large R) flux approaches the steady state result, 1/2. It is interesting 
to note that although the main change in the flux seems to take place over a wide region in R, the dynamical view 
with excitation localized to Rc works well, as demonstrated by us earlier in Ref. JT^ . 

In Fig. I we show a comparison between the MCWP results, the diabatic OBE results (D-OBE) and adiabatic 
OBE results (A-OBE). The agreement between the MCWP and A-OBE results is very satisfactory, whereas the 
D-OBE results fail by an order of magnitude for T = 0.3 mK. This failure increases further as T decreases, as 
shown in Refs. ]lO| , p^ . As discussed in the previous paragraph, the difference between D-OBE and A-OBE results 
suggest strongly that the basic condition for the validity of the local equilibrium model is not fulfilled for typical trap 
parameters. 

We compare the complex potential method and the Landau-Zener approaches to the MCWP results in Fig. ||. The 
LZD method saturates when the Landau-Zener excitation probability Plz becomes unity; until then all methods seem 
to agree well. However, beyond the saturation of the excitation the complex potential approach fails utterly. One 
should note that the saturation of the dynamical excitation is not the same as the saturation of the atomic excitation 
(also, the atomic excitation saturates to 1/2, but the dynamical excitation to unity). The complex potential method 
fails utterly when one approaches the saturation limit. The LZDD method agrees well with the MCWP results. 

We have used the A-OBE and LZDD method to calculate the flux for various detunings, and the results are given 
in Figs. H and 0. As can be expected, the saturation moves to larger fl when A increases. In Fig. 1(a) we start to see 
the signs of the failure of the Landau-Zener model at small T (and small A). 



VIII. CONCLUSIONS 



In this paper we have derived the adiabatic optical Bloch equations. When applied to the standard two-state model 
for loss of laser-cooled atoms from electromagnetic traps, these equations prove to be a fast and adequately accurate 
method to predict probabilities to reach any internuclear distance on the excited quasimolecule state. The A-OBE 
results match very well with the MCWP results, and also agree with the results from the qualitative LZDD model. 
The latter agreement is surprisingly good, but this can be attributed to a fortuitous definition of the edge of the 
interaction region. 

Although the A-OBE method is a useful tool, we are still far from having a proper general treatment of trap 
loss at small detunings. Firstly, the bound state structure of the attractive excited state is not taken into account. 
At small detunings the vibrational states associated with the attractive molecular potential overlap strongly and at 
strong fields they are also power broadened. In the language of time-dependent approach, the first "vibration" of 
the quasimolecule dominates over all the others. We can remedy the situation to some extent by applying the single 
passage result to multiple passages, as has been done e.g. in Ref. p6| . 

Secondly, the A-OBE method does not allow for momentum change; the motion on the ground state is given by 
the ground state velocity factor. When the excited state population of the quasimolecule decays back to the ground 
state, its kinetic energy distribution has been affected by the acceleration on the excited state. If this population is 
excited by e.g. another laser, this excitation depends on the new kinetic energy distribution. This effect is important 
in the cases of radiative heating jl^ and the radiative enhancement [ p7|j2^ ] . Currently only the MCWP method can 
handle the kinetic energy changes correctly [Tl| ]. 

Thirdly, the velocity factors diverge at classical turning points. Thus the A-OBE method is not capable of treating 
the case of optical shielding [p|JT^ , which involves excitation to a repulsive quasimolecule potential by a blue-detuned 
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FIG. 4. The excited state flux Je at R = 143 oo as a function of the laser-induced coupling f2, calculated using the MCWP, 
D-OBE and A-OBE methods. Here A = 7at, and the other parameters are as indicated in the figure. 
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FIG. 5. The excited state flux Je at R = 143 ao as a function of the laser-induced coupling Q, calculated using the MCWP, 
complex potential, LZD and LZDD methods. Here A = 7at, and the other parameters are as indicated in the figure. 
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FIG. 6. The excited state flux Je at i? = 143 oo as a function of tlie laser-induced coupling Q, calculated using the A-OBE, 
LZD and LZDD methods for T = 0.3 mK. 
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FIG. 7. The excited state flux Je at 7? = 143 oo as a function of tlie laser-induced coupling Q, calculated using the A-OBE, 
LZD and LZDD methods for T = 1.0 mK. 
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laser. This problem might be solved by replacing the WKB wave function by a proper Airy function solution in our 
derivation of the A-OBE method |p9|. 

Finally, we have given here as an example only the case of one ground state and one excited state. In reality we 
have several states involved in the cold collision process, e.g., the various partial waves and hyperfine states. In such a 
situation it is not so easy to write out the transformation to the adiabatic representation in an analytic form. However, 
there are various methods for doing the change of basis numerically. A problem may arise from the fact that we had 
to redefine the velocity factors using physical arguments — it is not obvious if such a redefinition in a multistate case 
would be as straightforward and clear as in the two-state case. In any case, the A-OBE method should be capable of 
handling adequately the case of several partly overlapping strong crossings in a system of close-lying states. In such a 
situation the Landau-Zener methods are expected to fail — the A-OBE method can help in testing the validity of the 
Landau-Zener methods in nontrivial crossing situations. 

For small detunings the Condon point moves to very large distances and the quasimolecule potential does not 
have the simple dependence any more. Furthermore, the retardation effects become important, and then the 

quasimolecular line widths become R dependent even in the diabatic representation. The A-OBE method should be 
very useful in studying these situations, as it is much faster than the MCWP method when exploring a wide range of 
cases with varying laser parameters (A and il) and quasimolecule potentials (C3, I). 

Despite some limitations the A-OBE method is a clear leap forward in treating theoretically the cold collisions in 
light fields. The discussion above, however, should be taken as a warning against trusting the method too blindly. 
The results given in this paper and in our previous report indicate nevertheless that the method is very good in 
predicting the behaviour of trap loss due to near- resonant light. 
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